Dispersion relations for circular single and double dusty plasma chains 
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We derive dispersion relations for a system of identical particles confined in a two-dimensional 
annular harmonic well and which interact through a Yukawa potential, e.g., a dusty plasma ring. 
When the particles are in a single chain (i.e., a one-dimensional ring) we find a longitudinal acoustic 
mode and a transverse optical mode which show approximate agreement with the dispersion relation 
for a straight configuration for large radii of the ring. When the radius decreases, the dispersion 
relations modify: there appears an anticrossing of the modes near the crossing point resulting in a 
frequency gap between the lower and upper branches of the modified dispersion relations. For the 
double chain (i.e., a two-dimensional zigzag configuration) the dispersion relation has four branches: 
longitudinal acoustic and optical and transverse acoustic and optical. 

PACS numbers: 52.27.Lw; 63.20.Dj; 64.60.aii 
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I. INTRODUCTION 

An interest in the structural and dynamical properties 
of ordered many-body systems confined by various ge- 
ometries is due to the recent advances in studies of, e.g., 
confined plasma [l|-13]- Knowledge of the properties, e.g., 
the ground state configurations and the normal mode 
spectra of the system allows one to calculate thermody- 
namic characteristics such as the free energy, entropy, 
etc. For example, theoretical investigation of properties 
of one- and quasi-one dimensional linear infinite chains 
of particles were presented in Refs. 0-01 • Experimen- 
tal studies of vibrational properties for one-dimensional 
linear finite chains of particles were performed in Refs. 

The ground state configurations and the normal mode 
spectra of a quasi-one-dimensional multichain plasmas 
of charged particles interacting through a screened 
Coulomb potential were investigated for monodispersed 
particles [3, 0] and for binary mixtures [7| . The one- to 
two-dimensional continuous zigzag transition was studied 
experimentally [2'| for small clusters of particles interact- 
ing through a Yukawa potential, and theoretically in [Sj] 
for confined along one direction systems of strongly in- 
teracting particles. Continuous and discontinuous struc- 
tural phase transitions in confined systems were also in- 
vestigated [y, 01 ■ The evidence of a direction-dependent 
melting and reentrant meltin g as a function of the density 
of particles was found in [y, llll | . In the experiment [9| 
a transverse optical mode in a one-dimensional Yukawa 
chain which had negative dispersion (i.e., when phase and 
group velocities were oppositely directed) was observed. 
The dynamic structure factor of one-dimensional chains 
with different types of disorder was studied in [12, [l3 ■ 
The linear and weak non-linear wave propagation in one- 
dimensional lattices of dust plasma was investigated in 
Ref. [1J|. The studies of the energy spectrum, the 
eigenmodes, the density of states and melting proper- 
ties of finite two-dimensional Yukawa lattices confined in 



parabolic and/or hard-wall external potential were per- 
formed in Refs. [ni[il-[i3. 

The properties of two-dimensional Coulomb clusters 
confined in a harmonic potential were theoretically stud- 
ied in Ref. 18]. A numerical analysis of normal modes 
of the ground-state configuration was performed, and the 
stability of the ring cluster was analyzed. In 3] one- 
dimensional and quasi-one-dimensional strongly-coupled 
dusty plasma rings were created experimentally, and the 
longitudinal and transverse dispersion relations for the 
one-dimensional single chain were measured. The mea- 
sured dispersion relations were in a good agreement with 
theoretically predicted results for one-dimensional ring 

In this work, we theoretically investigate the normal 
mode spectra for single and double circular chains. We 
analyze these spectra as a function of varying parameters 
such as the radius of the channel a and the density of 
particles n — N/2TTro, where rg is the equilibrium radius 
of the chain and N is the number of particles in the 
chain. In addition, we compare the obtained spectra for 
a double circular chain with those for a single chain (i.e., 
in the vicinity of the zig-zag transition) and for a double 
straight chain (i.e., for large radii of the circular chain). 

The paper is organized as following. First we discuss 
the normal mode spectra of a single circular chain in Sec. 
II. Then in Sec. Ill we derive the dispersion relations for 
an infinite double straight chain and finally, we present 
the normal mode theory for the case of double circular 
chain in Sec. IV. Conclusions are given in Sec. V. 



II. SINGLE CIRCULAR CHAIN 

Let us first consider a system of N particles in a cir- 
cular parabolic confinement which interact via a Yukawa 
potential. Fisrt, we discuss the case of a single chain. The 
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FIG. 1; The dispersion curves for single circular (triangles 
for longitudinal modes and squares for transverse modes) and 
straight (lines) chains for two different values of density n — 
0.596 (curves 1) and 0.892 (curves 2) (cp. Ref. 0) and for 
the circle radius a — 16. Dashed lines are guides for eye. The 
same symbols are used in Figs. 2 and 3. 



potential energy of the system in dimensionless form is: 



C/ = ^(.-n)VE' 



(1) 



i<j 



where a is the radius of the circular channel, r^ is the co- 
ordinate of ith particle, r.y is the interparticle distance, 
and K is the inverse Debye screening length. The trans- 
formation to dimensionless units is defined by the follow- 
ing relations: r^j = fij/ro, a.^j = a/ro, n = n/ro, k = 
kro, U = U/Uo, and finally U^ — aq^ je" , r\ = q^ /at, 
where all the dimensional units are labeled with tilde, 
and a is the strength of the circular parabolic confine- 
ment, q is the charge of a particle. The equilibrium chain 
configurations required for the calculation of the normal 



modes are found by minimizing the potential energy of 
the system. 

The Lagrangian of this system in polar coordinates is 



n 



(2) 



InC^ 



a,f3,n — n' 



where (a, /3) = (r, </>) are the coordinate indices, (n, n') — 
1..N/2 are the indices of the cells, g^ = r„ — rg, <?„, — 
4'n' — (f'n'o are the coordinates, and vq is the equilibrium 
radius of the chain. The coefficients of the elastic cou- 
pling are given by: 



X''f'{n-n')^dgo,d^,UUg, 



(3) 



where the eq denotes the set of the equilibrium values of 
all coordinates r„o, 0no- Solving the problem of normal 
modes, we obtain: 



A = 



J^cl>(j> ^01- 



r'n 



Ar^ \rr j ""^ I Q i 

where the matrix elements are defined by 



(4) 



(5) 



Here j — n — n' , Aip = 2tt/N , k ~ 0..7r/a is the wave 
vector, where a — tqA^ is the distance between nearest 
neighbours. 

The solution of the problem of normal modes in po- 
lar coordinates results in the following general dispersion 
relation for a circular single chain of particles: 



J 



1 I'e^^le— 



TT? ^ 2^2^ ~^ ^^^^^ + Kr + 1) 1^=^^ cos {jA(j)) (1 - cos (krojAcj))) 

3 = 1 







± 



[JV/2] g„«. 



1 - Y^ -^- (kV^ + Sur + 3) |r=r, (1 - cos (krojAcj))) 



i=i 



y^ — — (k^t^ -I- Kr -I- l) \r=rj sin {jA(p) sin {krojAcp) 



(6) 



sin (^) and Acj) is the angular dis- In the hmiting case of large radius ro ^- oo (Acj) -> 0) 



where rj = 2ro 

tance between nearest neighbours. 



Eq. ^ can be reduced to the analytical result for the 



single straight chain (cp. Ref. Q): 






{kv + 1) |r=r, (1 - COS (kja)) , (7) 



E' 



(^{nr + if + l) 1,=,^ (1 - cos (fcja)) . (8) 



where r, 



ja. 



Relation © can be further simplified if we consider 
nearest neighbours approximation (NNA). In this case, 
one should omit all the terms in the sums of Eq. ([5]) 
except for the term j — 1 which defines the elastic coef- 
ficient X between the nearest particles. Then we obtain: 



X cos (A(/)) (1 - cos (fcA(/>)) ± ^J{l - X (1 - cos {kA(j)))f + x^ sin^ (A0) sin^ (fcA^). 



(9) 



r 



The dispersion curves calculated using Eq. ^ are pre- 
sented in Fig. 1 for the same parameters (for densities 
n = 0.596; 0.892 and corresponding numbers of parti- 
cles N = 60; 90, the radius of a parabolic channel is 
a = 16) for single circular (open circles) and straight 
(lines) chains. Results were checked by computing nor- 
mal mode spectra using the method of Ref. [3j] and found 
to be in good agreement with the calculations presented 
here. It is worth noting that the agreement between the 
dispersion curves is fairly good not only for large chains 
but remarkably even for chains with rather small radii as 
it can be seen from Fig. 2. The small difference in the 
dispersion curves for a circular chain and for a straight 
chain is related to the fact that there is a non-zero cou- 
pling between the directions of oscillations in the trans- 
verse (radial) and the longitudinal (tangential) modes in 
case of a circular chain (Fig. 1 and Fig. 2). For a system 
with non-zero curvature these modes are not indepen- 
dent unlike for the case of a straight chain. In order to 
clarify the dependence of the dispersion curves on the 
curvature of the channel, we show in Fig. 2 these curves 
for different radii a but for the same density n — 0.84. 
As it can be seen from Fig. 2, the difference between the 
curves for circular and straight chains diminishes when 
the radius of the chain increases. The main qualitative 
difference is observed for intermediate A:- values while for 
small and large k the dispersion curves are approximately 
identical. Although we should also note that there is a 
shift of the longitudinal branch in the area of interme- 
diate and large k which increases with decreasing the 
channel radius a. In particular, we observe the appear- 
ance of a "gap" between the transverse and longitudinal 
branches of the dispersion curves ( "anticrossing" ) for a 
circular chain. Note that the modes are discrete (see, 
e.g.. Fig. 2(a)) and therefore this gap is well defined only 
for large enough density n when the dispersion curves 
become quasi-continuous. It is also seen from Fig. 2 that 
the width of the gap increases with decreasing radius of 
the circular chain tq. Inside the frequency gap no oscil- 
lations can be excited in the system. 

Now we will discuss the dispersion curves as a func- 
tion of the density of particles n for small channel radius 
cr = 4. The corresponding plots are shown in Fig. 3. 



Figs. 3(a) refers to low densities n which is equivalent to 
the situation of weak interparticle interaction. The dis- 
tance between the nearest particles is large enough, there- 
fore in spite of the strong coupling between transverse 
and longitudinal modes of the dispersion curves (shown 
by symbols in Fig. 3(a)) are similar to the curves for a 
straight single chain (shown by lines). For higher den- 
sities n (Figs. 3(b)-(c)) these differences become more 
pronounced. For high density n, particles are situated 
closer to each other and therefore they interact stronger. 
This results in an anticrossing of the transversal and lon- 
gitudinal dispersion branches (lines). Due to the circu- 
lar symmetry of the system, there is a coupling between 
transversal and longitudinal modes that leads in turn to 
the splitting of the transverse and longitudinal branches. 
This coupling is absent in case of a straight chain there- 
fore there is no gap in the normal mode frequencies (see 
Fig. 3 (lines)). It should also be noted that the transverse 
branch gradually approaches zero at the edge of the band. 
As soon as the longitudinal frequency approaches zero, a 
phase transition to the double-chain configuration takes 
place [Sj. The condition for the zig-zag phase transition 
was discussed in 0, Q . 



III. STRAIGHT DOUBLE CHAIN 

The dispersion properties of a quasi-one-dimensional 
classical Wigner crystal were investigated in detail in 
Ref. Q in the presence of dissipation. The results with- 
out dissipation were analyzed in Ref. [y] and correspond- 
ing result for the binary system of particles in [fl\ . Never- 
theless it is worth to revisit here the dispersion relations 
for a multichain case. 

Let us consider a system of particles confined to a 
straight channel. The potential energy of the system in 
dimensionless form is described by: 



u-Ev^ + E- 



(10) 



i<j 



Solving the problem of normal modes for this system re- 
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FIG. 2: The dispersion curves for single circular (open circles) 
and straight (lines) chains for density n — 0.84, for various 
values of the radius a = 2.94 (a), 6.0 (b) and 12.09 (c). 
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FIG. 3: The dispersion relations for circular (open circles) 
and straight (lines) chains for a channel of radius a — 5 and 
for different values of the density n = 0.57 (a), 0.69 (b) and 
0.88 (c). 



suits in the following dynamical matrix: 



A 



where the matrix elements are given by 



A XX 
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Aif\ 




/I 
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\xx 
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\yy 





\yy 1 
^^22 / 




\0 1 



Kf = lT.^tfU)' 



^ikja 



(11) 



(12) 



where j denotes the jth elementary cell of the chain, 
(I, s) — (1, 2) is the index of the chain, a is the interpar- 
ticle distance (in the same chain) and (a, /3) — (x, y) are 
the coordinates. 

The independent matrix elements of the dynamical 
matrix A are given in Appendix A. Other elements of 



the dynamical matrix (llip can be defined using the fol- 
lowing relation from the theory of normal modes [19j : 



A = At. 



(13) 



It is clear that the presented results reveal some differ- 
ences with similar results obtained in Jsi, 6] for a straight 
zigzag configuration. Let us briefiy discuss them. Co- 
efficient A^ defines the elastic coupling between the a 
coordinate of some particle of the Zth chain from zeroth 
cell and /3-coordinate of a particle of s-th chain from 
jth cell (see Fig. 4). (Further we will call these coeffi- 
cients "a/3-couplings".) As it follows from Refs. [a, Q all 
a/3 = xy(yx) couplings are equal to zero. Our calcula- 
tions give the same result only for particles of the same 
chain, i.e., when I = s. However, in the general case 
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FIG. 4: Scheme illustrating elastic couplings in a straight 
double chain. Coefficients Af^ describe coupling between ^- 
oscillations of Zth chain and ^-oscillations of sth chain. Plots 
(a) and (b) show the couplings between longitudinal (a) and 
between transverse (b) of the same or different chains. Plots 
(c) and (d) show the couplings between the longitudinal and 
transverse oscillations in different chains [see Eq. (Ill|l ]. 



I 7^ s these matrix elements are non-zero. This can be 
demonstrated by a direct calculation and be explained 
from physical arguments. 

First, we discuss the case of a/3-couplings for / — s 
(see Figs. 4(a)-(b) for Af^ and Xfl). If some oscilla- 
tions are excited in the y-direction in the system then 
this would obviously result in a change of the length of 
couplings (and as a result, in its strength) with other 
particles and, as a consequence, this induces oscillations 
of other particles in the x-direction (see Fig. 4(a)-(b)). 
However, within the harmonic approximation this varia- 
tion of coupling strength is of second order in the strength 
of coupling or first order in the y-deviation. Thus we can 
conclude that the oscillations cannot be excited in this 
case as verified by direct calculations. The application of 
the same procedure to the yx-coupling leads to a similar 
result. 

However the situation is different for the case I ^ s 
(see Figs. 4(c)-(d) for A^| and A^j)- In this case, varia- 
tions of the strength of the a/3-coupling between particles 
from different chains (which take place for, e.g., trans- 
verse oscillations) are of the same order as the deviation 
from equilibrium. A displacement along the y-direction 
of a particle results in the corresponding displacement 
of the nearest particles of the adjacent chain in the x- 
direction (and vice versa). It is reflected by the fact 



that the corresponding coefficients (A^f and A^2 ) ^^^ non- 
zero, although the absolute values of these coefficients are 
smaller than those for the xx— and yy-couplings. 

Therefore the results presented in Refs. [a, Q sxe valid 
only for a limiting case of a very narrow double chain (i.e., 
when d << a). In this case the a/3-couplings for I = s be- 
come approximately applicable to the case I y^ s. This re- 
mark remains valid for all other multichain configurations 
considered in Refs. |5|, |6|. Thus presented formulas 
(IA1[) - (IA5|) converge to the ones obtained in Refs. 
when the condition - <C 1 is fulfilled. In this case formu- 
las dll]), (|Ai|)-(|A4| are valid and limd«a A"f\a^p,i^s = 0. 

Let us now discuss the dispersion relations for a double 




FIG. 5: The dispersion curves for a double straight chain 
for various values of density n and the equilibrium interchain 
distance d: n ^ 0.90; d = 0.209; a = 4.5; iV = 26 (a)-(b) and 
n = 1.19; d = 0.748; a = 5.0; N = 38 (c)-(d). The dependen- 
cies were calculated using relations (jllI) . (IA1|) - (|A5|) [(a) and 
(c)] cases and using the same relations but with zero second 
diagonal in the dynamical matrix ^] [(b) and (d)]. 



straight chain presented in Fig. 5. As shown in Figs. 5(a)- 
(b) for smaU ratio | = 0.081 (- < 1) the curves cal- 
culated using Eqs. (fTTI) . ()A1|) - ()A5|) and formulas from 
Ref. [6( look similar. Fig. 5(b) shows that the LA-branch 
(longitudinal acoustic) of vibrations intersects both the 
transverse branches (TA and TO). This fact is a conse- 
quence of neglecting the interaction between the trans- 
verse and longitudinal vibrations in the dynamic matrix. 
Taking into account the interaction between transverse 
and longitudinal vibrations leads to the "anticrossing" in 
the crossing points of LA and TO, LA and TA branches 
as shown in Fig. 5(a). As a result, the LA, TA and TO 
branches change and acquire frequency gaps. On the 
other hand, the difference between these two cases (i.e., 
with and without the xy— coupling taken into account) 
becomes more pronounced when the ratio - is not small: 
f = 0.389 (Figs. 5(c)-(d)). From Fig. 5(c) it is seen that 
only the LO and TA branches experience hardening (un- 
like in Fig. 5(d) where all the branches except the LA 
branch show only a weak dependence on the wave vec- 
tor). Fig. 5(c) shows that for density n = 1.19 there 
is a narrow frequency gap between the longitudinal and 
transverse branches. 



Thus the account of the interaction between the trans- 
verse and longitudinal oscillations for different chains 
in case of a circular double chain is important and 
leads to an "anticrossing" in the corresponding disper- 
sion branches (Figs. 5(a), (c)). 




FIG. 6: The dispersion curves for double circular (the left 
column) and straight (the right column) chains for various 
values of channel radius but for the same density n — 0.934 
close to phase transition point. The channel radius and the 
number of particles is cr = 5.26, A'^ = 30 - (a) and a — 26.66, 
iV = 150 - (c). 




IV. CIRCULAR DOUBLE CHAIN 
A. Dispersion relations 

Next we study the dispersion properties of a double 
circular chain. The potential energy of the system is 
given by ([Ij. Due to the symmetry of this system the 
Lagrangian can be written in polar coordinates in a fol- 
lowing manner }19| : 



l.n 



L- ?E (cf + M. 



(14) 



1 Y: Kfin^n')q:„C. 



where {I, s) = (1, 2) are the indices of the chains, (a, /3) = 
(r, (f>) are the coordinate indices, gj„ — Vsn — fgo, qf„, = 
4>in' — 4>in'o are the coordinates, r;o is the equilibrium ra- 
dius of the ^th chain, and {n, n') = 1..N/2 are the indices 
of the cells. The coefficients of the elastic coupling are 
given by: 



A"f (n - n') = dq^^d U\eq, 



(15) 



where eq denotes the set of coordinates rsno , 4'sno at their 
equihbrium value. Solving the problem of normal modes, 
we obtain: 
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where the matrix elements are defined by 



Ar/-^E^r/o-)e^"^°^'^'- 



(16) 



(17) 



FIG. 7: The dispersion curves for double circular chain for 
various values of the channel radius a = 11.33 (a), (b); 11.30 
(c), (d) and 11.20 (e), (f) near the single-to-double ring tran- 
sition point. Comparison of the dispersion curves obtained 
analytically using Eq. H16|l (right column) and those calcu- 
lated numerically (left column). The number of particles for 
all cases is the same: A^ = 64, and the corresponding densities 
are n = 0.896 (a), (b); 0.898 (c), (d) and 0.906 (e), (f). 



Here j — n — n' , and I, s are defined above, N is the 
number of particles, A(f) — Air /N and k is the wave vec- 
tor. The independent matrix elements of the dynamical 
matrix A are given below in Appendix B. Taking into 
attention the asymptotic relations given in Appendix B 
in the limit r;o -^ oo, (A0 — > 0) and - ^- we obtain, 
as a consequence, the following relations: 



limA-^Af, 



limAr,^ = Aff^O, 



hm A- = Af/ = 0, ? ^ s hm A[^^ = Af; = 0, / ^ 



lim ^A"^"^ 
wo 



A XX 
II 1 






A XX 



0, I ^ s lim ^-A 

^10 



0r 



Af/ 



0, l^s. 



(18) 



By direct calculation one can show that the obtained 
formulas (IT6l) . (|B1[) - (|B7I) converge to the ones found in 
the previous section, (fTTj) . (|A1|) - (|A5|) [in the limit de- 
fined by expressions (|B8p and p^ ]. and in the limiting 
case d << a to the corresponding results of 0. This is 
also verified by a comparison of the dispersion curves for 
double circular (left column) and straight (right column) 
chains which are presented in Fig. 6. The dispersion 
curves were calculated for various values of channel radius 
a but for the same density n sa 0.934 close to phase tran- 
sition point. The channel radius and the number of parti- 
cles are a = 5.26, A^ = 30 for (a) and a = 26.66, N = 150 
for (c) plots. Comparison of Fig. 6(a) and (b) shows the 
presence of a gap in the dispersion relations for a circu- 
lar chain and a straight chain. In addition, as we can 
see from Fig. 6(a), the LA-branch for small fc (fc w 0) is 
characterized by w sa which means that in the vicinity 
to fc = the longitudinal mode is pure rotational and it 
requires zero energy to excite. Comparison of Fig. 6(c) 
and (d) reveals significant similarity between the disper- 
sion curves for double circular and straight chains. This 
is reasonable because this case corresponds to small cur- 
vature. Thus we can conclude that differences between 
the dispersion relations of the circular and straight chains 
increase with decreasing channel radius a. 



B. The single- to double-chain transition 




FIG. 8: The dispersion curves for double circular (left col- 
umn) and straight (right column) chains for various values of 
density n beyond phase transition point (0.94(N=30)-(a),(b); 
1.13(N=36)-(c),(d); 1.25(N=40)-(e),(f)) for channel radius 
a = 5.0. 



Here we analyze the special case of the single-to-double 
chain transition for a circular chain. For that purpose, 
we considered the limit d — >■ in the expressions for a 
double circular chain. On the other hand, in order to 
verify the obtained results for this limit, we performed an 
independent numerical calculation for a double circular 
chain just above the single-to-double transition. The re- 
sults obtained using both methods are shown in Fig. 7 for 
a ~ 11.3 and n « 0.9. Comparison of numerical and an- 
alytical calculations reveals an excellent agreement (see 
Fig. 7). 

Note that at the zig-zag transition point, the length of 
the unit cell increases from a to approximately 2a (for 
each chain in the double-ring configuration). Therefore, 
the Brillouin zone shrinks (left column) giving rise to the 
new brunches in the dispersion relations. 



C. Role of the density 



fore, the best agreement between the dispersion curves 
for double circular and straight chains is reached for the 
lowest value of particle density n just above the critical 
density at the phase transition point. However this fact is 
not obvious: it is clear that the best correspondence be- 
tween the dispersion curves for double circular and dou- 
ble straight chains is reached when the unit cells for both 
chains are very similar. On the other hand, increase of 
the density n results in enhancing this similarity (due 
to decreasing angular spacing A0) and simultaneously in 
its diminishing (due to increasing interchain spacing d). 
Therefore we studied the imbalance between the "lattice 
constants" in the inner and outer chains, a2 — ai — d/S.<f> 
as a function of the density n and we found that the max- 
imum similarity in the unit cells for double circular and 
straight chains is indeed reached for critical density ric. 



The dispersion relations are shown in Fig. 8 for the 
same channel radius a — 5.0 for double circular channel 
(left column) calculated using Eqs. ([16]), ((BT|) - (IB7)) and 
straight (right column) calculated using Eqs. pTj) . (lAip - 
(|A5|) for densities n slightly above the critical density 
Tie at the single-to-double chain transition. As one can 
see from Fig. 8, difference between the dispersion curves 
for double circular and straight chains becomes more 
pronounced with increasing particle density n. There- 



From Fig. 8(a), (c) and (e) it can be seen that the 
gap between dispersion curves first increases with grow- 
ing density and then decreases. In addition the frequency 
of TO branch experiences hardening for small fc « as 
it is seen from comparison of Fig. 8(a) and (c). 

Thus Fig. 8 reveals that there is essential coupling be- 
tween LA and TO and LO and TA branches of the eigen- 
modes. 



V. CONCLUSIONS 

We derived an analytical expression for the dispersion 
relations for a single circular chain of interacting parti- 
cles and analyzed the obtained results for different radii 
of the chain and the influence of the curvature of the 
confinement channel. For large radius of the chain our 
results converge to the dispersion relations of a straight 
single chain [5|. 

An improved expression for the dynamical matrix for 
the case of a straight double chain was obtained that 
takes into account all the elastic (harmonic) couplings 
between the oscillations of the particles in the different 
chains. Taking into account the interaction between the 
transverse and longitudinal oscillations of the different 
chains for the case of a double straight chain leads to a 
weak transverse-longitudinal "splitting" (i.e., "anticross- 



ing"). 

We found an analytical expression for the dispersion 
relations of the normal-mode oscillations for the case of 
a double circular chain of interacting particles, and we 
compared these relations with those for a double straight 
chain Q . We showed that there are small differences be- 
tween the dispersion relations of a circular and a straight 
both a single and a double chains. These differences grow 
with increasing curvature of the parabolic confinement. 
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Appendix A 

Here we present the analytical expressions for matrix elements of the dynamical matrix A for a straight double 
chain (jlip discussed in Sec. III. The corresponding independent matrix elements of the dynamical matrix A are: 



A XX 
a 



E- 



{^{nr + lf 



a^ - {kt + 1) S 



|r=r_ H 3- ( (Kr 4^ 1) + 1 j \r=ja (1 " COS [kjo)) 

(Al) 



Af,^ 



1+E 



j=i 



[[kt 



\\d' 



{nr + 1) ( i - - 



(nr + 1) \r=ja (1 - cos (kja)) (A2) 



A XX 
Is — 



E' 



((Kr + l)Vl) h-^J a^-{Kr + l)d^ \r=r^cos( k(j 



l^s 



(A3) 



^Vs 



sr^ e 



j=i 



((Kr+l)' + l)d2-(«;r + I) b'-ij a^ |r=r„ cos ^ j - i j 



I ^ s (A4) 



p~ r / 
Af/ = Af; = z 5] ^ [[{nr + if + [nr + l) + l 



r=r_d\j - -\asinikij - -]a] , l^s (A5) 



where r_ — \ [j — ■k)^a'^ + d"^ and d is the equilibrium distance between the chains. Other matrix elements can be 
calculated as hermitian conjugates using Eq. (|13p . 
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In this Appendix we present analytical expressions for matrix elements of the dynamical matrix A (|16p discussed in 
Sec. IV for a circular double chain. The independent matrix elements of the dynamical matrix A for circular double 
chain are: 



AS* = i: 9^ [((-"• + D' + 1) »■>' ((i - A*) 

+ 5] r; 



J - 2 ) ^-^ 



iO ^5 



{nr + 1)" + 1 sin^ (jAcf)) + (nr + 1) (1 - cos ijA(j))y |^=-„ 



X (1 — cos (krojAcI))) 






'-+>''+>i"i^+^- 0-^'^* 



1 + cos^ 
r 



.? - 2 I ^'^ 



+ (Kr + l)sin^ ( [j--]A(l) 



10 ^5 



((Kr + 1)' + l) (1 - cos {jA(j))f + {nr + 1) sin^ (jA0)] |,=^,„ 



X (1 - cos (krojAcj))) 



(Bl) 



(B2) 



Af/ = --V^^oE^f ((- + !)' + !) --'((.?■ 



+ (Kr + 1) 1 I cos 

\\r2Q no 



A0 



1 + cos-^ j 



1 



A(P 



A(t) 



1 



=jr„ cos Zero i - - A(/) , l^ s 



(B3) 






-?'10''20 E 



J = l 



1 + cos^ [[j--]Acj, 



(«;r + l)sin^ j - - A(/) 



<" + "' + 'l"^ + S)"<(-5'^* 



|r=ro cos ( fcro ( i - - ) A(/) ) , / 7^ S 



(B4) 



p r 

Ah^ - «rfo E ^:^ (1 " ^°^ (•?'^'^)) ^i^ (•?'^'^) [('^'^ + 1)' + ('«'' + 1) + 1 lr=no sin (A:rojA0) 



(B5) 



Ai^ = irao E sin ( ( i - X ) A(?!. 



i=i 



{kt + 1) r2o ( r2o - rio cos ( ( j - - ) A^ 



+ ( (Kr + 1)^ + 1] rio ( rio - r2o cos ( ( j - - ) A0 



|,.=Fo sin ( Zero ( j - - ) A</) ) , ? 7^ s 



(B6) 



10 



^12 = -^^W X! ^^'^ ( ( -'^ 9 ) ^'^ 



J = l 



(Kr + 1) rio ( rio - r2o cos ( [ j - - ] Acj) 



(B7) 



(^(Kr + 1)^ + 1 j rso f 7-20 - r-io cos N j - - j A0 



ir=Fo sin ( Zero ( j - - ) A./) ) , l^s 



where r;o = v^r/ov^l — cos {jAcj}) and ro — y'^io + ''io ^ 2rior2o cos ((j — i) A(/)). Other elements of the dynamical 
matrix can be found in a similar way using Eq. (1131) . 

The asymptotic relations in the limit t/q -> cx), (A0 — > 0) and - — > 0, 




(B8) 



V2rlQ^/T^~cos{jA(i)) -^ ja, \rio - r^o cos ( ( j - - ) A0 ) | -> d, I ^ 



^^°^^°".^ + S)^°K0'4)^^ 



- 1 + cos^ [[j--]A(f> 



rior2o sm^ ( { j - ^ ] Ac/) 



J-2) « 
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